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Abstract: Five-dimensional scenarios with infinitesimally thin branes replaced by ap- 
propriate configurations of a scalar field were considered. A possibility of periodic extra 
dimension was discussed in the presence on non-minimal scalar-gravity coupling and a gen- 
eralized Gibbons-Kallosh-Linde sum rule was found. In order to avoid constraints imposed 
by periodicity, a non-compact spacial extra dimension was introduced. A five dimensional 
model with warped geometry and two thin branes mimicked by a scalar profile was con- 
structed and discussed. In the thin brane limit the model corresponds to a set-up with 
two positive-tension branes. The presence of two branes allows to address the issue of the 
hierarchy problem which could be solved by the standard warping of the four dimensional 
metric. Stability of the background solution was discussed and verified in the presence of 
the most general perturbations of the metric and the scalar field. 
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1 Introduction 

The idea of extra-dimensions of space-time inspired by the string theory received a lot 
of attention since the last decade or so as a possible solution to the problem of hierarchy 
between the Planck scale Mp\ ~ 10 18 GeV and the electroweak scale rriEW ~ 10 3 GeV [1-3]. 
The most attractive incarnation of the idea, proposed by Randall and Sundrum [3], involves 
one extra-dimension with non-trivial warped factor appearing due to the anti-de Sitter 
(AdS) geometry along the fifth-dimension. This is so-called RSI, the Randall-Sundrum 
model with two D3 branes on the S1/Z2 orbifold along the extra-dimension. In this model 
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the hierarchy problem can be addressed without introducing the large compactified volume 
of the extra-dimensions as suggested by Arkani-Hamed, Dimopoulos, Dvali (ADD) [1, 2], 
by the virtue of the non-trivial warped geometry along the extra-dimension. The model 
requires the presence of two singular D3 branes of opposite tension. There were many 
attempts to avoid the presence of thin (singular) branes. It has been shown [4-12], that the 
positive tension brane could be smoothed (or, in other words, regularized) by a background 
scalar field configuration which we term here as the smooth or thick brane. However there 
is no satisfactory simple strategy to model the negative tension brane at least for scalar field 
minimally coupled to gravity. For example of an existing attempts to generate negative 
tension branes see [13]. Also, as shown by Gibbons et al. [14], periodicity of set-ups 
like RSI is generically in conflict with the idea of a smooth non-trivial scalar profile. So, 
concluding, smoothing the RSI scenario is severely limited by periodicity and impossibility 
of generating a negative tension brane by a scalar field configuration. The conclusion holds 
at least for the case when the scalar field is minimally coupled to gravity, the generalization 
will be considered below. 

It has been pointed out by Randall and Sundrum in their second seminal paper that 
the extra-dimension can be infinite and yet it can lead to nearly standard 4D gravity [15]. 
The main idea in this second paper (RS2) is that a single D3 brane of positive tension 
is embedded in a 5D AdS geometry and the gravity is effectively 4D at large distances. 
The 4D graviton is localized on the brane which implies finite 4D Planck mass. Since a 
positive tension brane could be mimicked by a scalar field, therefore the RS2 seems to be 
more attractive from that perspective. However, the RS2 model having just one D3 brane 
is not addressing the hierarchy problem. Since solving the hierarchy problem is the main 
motivation for our study, therefore we will try to improve the scenario by introducing a 
second brane. 

The purpose of present work is twofold: 

• to see if one can overcome the above mentioned obstacles (periodicity and positivity 
of brane tension) to achieve a smooth version of RSI in modified gravity with the 
scalar field non-minimally coupled to the Ricci scalar, and 

• to verify if one can address the hierarchy problem with two thick branes (which in a 
certain "brane limit" mimic two positive tension singular branes) with non-compact 
warped extra-dimension. 

As far as the first question is concerned, we will show, by generalizing the Gibbons-Kallosh- 
Linde sum rules [14], that it is not possible to achieve the periodicity for scalar field and 
metric solutions even with the scalar field non-minimally coupled to the Ricci scalar. The 
consistency conditions have been discussed following the strategy of Gibbons et al. [14] 
in the modified gravity set-up [16], however the authors consider the scalar field in the 
bulk with singular branes, this is exactly what we want to avoid for our set-up, i.e., we 
wish to have smooth branes instead of singular branes. Another attempt to overcome the 
problem of periodicity was discussed in [17]. Concerning the issue of the positivity of a 
brane tension generated by a scalar profile we also find that even with non-minimal scalar 
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couplings there is no way to generate a negative tension brane. Therefore we turn our 
attention to models with only positive tension branes (e.g. the RS2) and non-compact 
(to avoid restrictions imposed by periodicity a'la [14]). Since we find a satisfactory model 
with a scalar that is minimally coupled to gravity, we restrict ourself to that scenario. 
In order to be able to address the hierarchy problem we will propose a model with two 
thick (smooth) branes, which, in an appropriately defined limit (so called brane limit) 
approaches two singular branes. The limiting version of the model was discussed earlier 
by Lykken and Randall (LR) in [18]. As we will show, in our set-up of two thick branes, 
different possible solutions for the warped factor can emerge, for instance, we can have 
the AdS or Minkowski geometry in different regions along the extra-dimension. We will 
discuss three such configurations, (i) the two thick branes between the AdS vacua so that 
we have warped geometry and hierarchy problem could be addressed in this set-up (this is 
the thick brane version of the Lykken- Randall model [18]), (ii) the case when we can have 
the Minkowski background in between the two branes and the AdS geometry to the right 
and left of both branes and, (Hi) with the Minkowski geometry to right or left of both 
branes and the AdS in the other regions along the extra-dimension, which gives the thick 
brane version of the Gregory, Rubakov, Sibiryakov (GRS) model [19] except that we have 
the second brane also with positive tension instead of the negative tension as in the original 
GRS model with singular branes. In all the above mentioned three cases we study stability 
of the background solutions and also existence and localization of the zero-modes of the 
scalar, vector and tensor (SVT) perturbations of the solutions. The issue of stability and 
localization of the zero-modes have been extensively discussed in the past in the context 
of thin as well as thick brane scenarios [7, 20-30]. We find that the SVT perturbation 
equations could be transformed into a supersymmetric quantum mechanics form so that 
they guarantee stability of these perturbations in all the configurations considered above 
and also the absence of tachyonic modes. It turns out that the zero-mode of the tensor 
perturbation wave function (that corresponds to the 4D gravitons) is localized in the cases 
(i) and (ii) but it is quasi-localized in the case (Hi) mentioned above. The zero- modes 
corresponding to the scalar and vector perturbations are not localized, as they are not 
normalizable modes, consequently they do not affect the 4D physics. 

The paper is organized as follows. Possibilities and difficulties of constructing smooth 
generalizations of the RSI in the modified gravity scenario is considered in Sec 2. In 
Sec. 3 the background solutions of the two thick brane set-up is discussed in the stan- 
dard Einstein-Hilbert gravity. Possible different geometric configurations for the two thick 
branes are discussed in Sec. 4 along with their brane limit solutions. Sec. 5 is dedicated 
to the linearized Einstein equations. The SVT decomposition of the perturbations and the 
gauge choice is discussed in Sec. 6. In Sec. 7 we have derived the equations of motion 
for the SVT perturbations. The issues of stability and localization of zero-modes of the 
SVT perturbations are covered in Sec. 8. The last section, Sec. 9, contains our conclu- 
sions. Appendices are dedicated to conventions and some mathematical details for the 
SVT decomposition of metric perturbations. 
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2 Thick brane generalization of RSI in modified gravity 



Our first goal is to mimic (regularize) D3-branes which appear in various five-dimensional 
(5D) scenarios that solve the hierarchy problem by warping the metric along the extra 
dimension in the spirit of [3]. The most natural approach is to introduce a 5D scalar 
field (ft with a non-trivial profile (that satisfies equations of motion) that in certain limit 
could mimic a brane by approaching a delta-like energy distribution along extra dimension. 
However, as it was shown in [14], in the case of compact extra dimensions the idea of a non- 
trivial scalar profile (a thick brane) is severely restricted by the requirement of periodicity. 
Arguments adopted in [14] apply for a scalar that is minimally coupled to gravity. Therefore 
here, we are going to discuss first a class of models allowing for non-minimal scalar-gravity 
coupling: 

S = J dx 5 ^ lf(cft)R - \g MN V M <ftV N <ft - V{<P)} , (2.1) 

where f(4>) is a general smooth positive definite function of the scalar field (ft 1 , which 
is supposed to compose the D3-branes that are present in the RSI scenario [3]. In our 
convention capital roman indices will refere to 5D objects, i.e., M,N,--- = 0,1,2,3,5 
whereas, the Greek indices label four-dimensional (4D) objects, i.e., /x, v, ■ ■ ■ = 0, 1, 2, 3. In 
the Eq. (2.1) Vm is 5D covariant derivative, for conventions and details see the Appendix 
A . 

In other words, branes would be made of the scalar field while other fields could be 
dynamically localized in certain regions of the 5D space, see for instance [31-34]. Thick 
branes in the presence of non-minimally coupled scalar was discussed earlier by [35, 36]. 

We will look for a solution of the Einstein equations with the following 5D metric 2 
ansatz 

ds 2 = e 2A( -^r]^dxV + dy 2 , (2.2) 

where r]uv(%) is the 4D metric and the warp function A(y) is only a function of the extra- 
spatial coordinate y. The Einstein's equations of motion, resulting from the action (2.1) 
are 

Run - \qmnR = \\ t mn + V M V N /(</>) - 5a/vV 2 /(^)| , (2.3) 

with the energy-momentum tensor for the scalar field eft as 

Tmn = V M cj)V N <ft - g M N Q(V0) 2 + V(<ftyj . (2.5) 



1 From the action (2.1) one can infer that (f>(y) and V((f>) have the mass dimensions 3/2 and 5, respectively 
and also we assume that the scalar field 4> has only dependence on the extra-spacial coordinate y. 
2 We choose the metric signature ( — h + + +). 
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From the Einstein equations (2.3) and (2.4), one can get the equations of motion for the 
metric ansatz (2.2) as, 




(2.6) 



(2.7) 



(2.8) 



where prime ' denoted the derivative with respect to y, and it is understood that / and V 
are functions of the scalar field 4>(y). 

2.1 Thick branes with periodic extra dimensions 

Here we would like to verify if the scenario with a non-trivial profile of a bulk scalar could 
be consistent with periodicity in the case of compact extra dimension. The authors of [14] 
derived elegant, simple and powerful sum rules that severely restrict thick brane scenarios 
with periodic extra dimensions. From our perspective the most relevant result obtained 
there is the following condition that must be satisfied for periodic extra dimensions with a 
bulk scalar <p when singular branes are absent: 



The above result implies that non-trivial scalar profiles are inconsistent with periodicity, 
the only allowed configuration is <p =const.. The sum rule (2.9) was obtained assuming 
minimal scalar-gravity coupling. In the following we are going to generalize the result for 
the case of non-minimal coupling described by the action (2.1). 

It is easy (subtracting Eqs. (2.6) and (2.7)) to derive an equation of motion that 
contains only the warp function A{y) and the input profile 4>(y). 



It is useful to rewrite the above equation by the change of variables X{y) = A'{y): 




(2.9) 



3fA" = f'A> - f" - ^'f. 



(2.10) 



X'(y) 



F(y)X(y) + G(y), 



(2.11) 



where, 



F{y) 



f'(y) 



(2.12) 



G(y) 



3/(v)' 

f"(y) + \4>'{y? 
3/(y) 



(2.13) 



We assume that the profile is periodic with a period L: 



4>{y + L) = 4>(y) 
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then f((f>) and consequently, F{y) and G{y) are also periodic with the same period L. Since 



£dyF(y)=0 



(2.14) 



it is straightforward to notice that the solution of the homogeneous part of Eq. (2.11) 



VI \ v f y F(s)ds 

is periodic as well. 

The inhomogeneous equation (2.11) could be rewritten in the following form 



(2.15) 



[Z(y)X(y)]' = Z(y)G(y), (2.16) 
where Z(y) = Z(0)e~ F ( s ) ds [ s a solution of the following homogeneous equation, 

Z'(y) = -F(y)Z(y). (2.17) 
Integrating (2.16) over the period we obtain the following condition: 



jdy G(y)Z(y) = 0, 



(2.18) 



which constitutes the proper generalization of the Gibbons-Kallosh-Linde sum rule (2.9). 
For F(y) defined in (2.12) one obtains Z(y) explicitly 



Z{y) = Z(0) 



1 1/3 



(2.19) 



Then, after integrating by parts, the sum rule (2.18) reads 

2 -i / // \ 2' 



dy 



4 1 



+ 



1 



3 /1/3 \ f J 2 V/2/3 



(2.20) 



The above sum rule again implies that even in the presence of non-minimal couplings, 
f(4>)R, only the trivial profile, <f> =const., is consistent with periodicity. 

Note that (2.20) holds also for multicomponent scalar fields, therefore even in that 
case non-trivial profiles in the absence of singular branes are excluded by periodicity. 

2.2 Negative tension brane in modified gravity 

In the RSI set-up one of the two 3-branes must have a negative tension. Therefore in this 
subsection we turn to the question weather a negative tension brane can be constructed out 
of a real scalar field in the modified gravity scenario when the scalar field is non-minimally 
coupled to gravity as in (2.1). By the virtue of the result of the previous subsection we 
discard the possibility of periodic extra dimensions. Let us assume, without loosing any 
generality, that the scalar filed 4>{y) has a kink-like profile, 



Hv) = ^tanh(/3y), 



(2.21) 
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where f3 is a brane-thickness controlling parameter. For f3 — > oo (the brane limit), as it will 
be discussed in details in Sec. 4, the profile 4>(y) generates singular energy density localized 
at y = that could mimic a D3-brane. 

The action for the kink configuration (2.21) can be written as, 



I 



(2.22) 



while the action for a brane localized at y = with a negative tension (A > 0) reads 



Sd3 



d 5 x^X5(y). 



(2.23) 



Using the equations of motion (2.6) and (2.7) we can rewrite the action (2.22) as follows 



dy 



+ V 



J dy l-^'f + \2{A'ff + 8A'f] • (2.24) 



As it will be clear from the next section, the only interesting set-up is such that the warp 
function A(y) reaches its maximum at the brane location (so A'(0) = 0), therefore among 
the above terms only the very first one contributes to the brane tension. However, as it 
is seen from (2.23) there is no possibility to reproduce the sign required by the negative 
tension. Therefore we conclude that a single kink-like profile can generate only a positive 
tension brane even in the case of modified gravity. 



2.3 Conclusions on thick brane generalization of the RSI model 

As we have shown in the proceeding subsections there is a conflict between the RSI scenario 
and the idea of branes generated by bulk scalar profiles: 

• As shown in Sec.2.1, even in the presence of the non-minimal scalar-gravity coupling 
f(cf))R, periodicity in the extra coordinate can not be reconciled with a non-trivial 
profile. 

• One of the branes in the RSI scenario must have negative tension, however as we 
have shown in Sec. 2. 2 even if scalars interact non-minimally with gravity there is no 
way to generate a brane with negative tension. 

The above observations prompt to give up compactness and therefore to discuss a possibility 
of mimicking the RS2 model with non-compact extra dimension. Since we would like to 
allow for the solution of the hierarchy problem by the virtue of warping the metric along 
extra dimensions, we will introduce a scalar field, the profile of which could mimic a scenario 
with two branes of positive tension with AdS metric in between them. This is what we are 
going to discuss in the following sections limiting ourselves to the case of minimal scalar- 
gravity coupling, however the analysis could be easily extended to non-minimal scenarios 
as well. 
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3 Two thick branes and the background solutions 

From now onwards we will adopt the following action for a 5D scalar field minimally coupled 
to the Einstein-Hilbert gravity 



S = f dx 5 ^9 |2M 3 i? - ^ mjv V m 0Vjv0 - V{<f>)\ , 



(3.1) 



where gMN is the warped 5D metric, 



9MN = I 1/' ^ ' 

The Einstein's equations and the equation of motion for <fi, resulting from the action (3.1) 
are 

Rmn - -^qmnR = 4jp TMN ' 

W-f = 0, (3.4) 

where V 2 is 5D covariant d'Alambertion operator (see Appendix-A for conventions and 
details) and the energy-momentum tensor Tmn for the scalar field <p{y) is, 



Tmn = Vm^Vat^ — gMN 



i(V0) 2 + y(<A) 



(3.5) 



From the Einstein equations (3.3) and (3.4), one can get the following equations of motion 
for the metric ansatz (3.2) 

24M 3 (A') 2 = ±(<t>') 2 -V(<t>), (3-6) 
12M 3 A" + 24Af 3 (A') 2 = -^(<P') 2 - V (</>), (3.7) 

0" + 4Ay-^ = O. (3.8) 

d(p 

We assume that the scalar potential V{4>) could be expressed in terms of the superpotential 
W(<j>) as [5], 

where W(cp) satisfies the following relations, 

and A'~^<«. ,3.0) 

Although the use of this method is motivated by supersymmetry, no supersymmetry is 
involved in our set-up. The method is elegant and very efficient, in particular it implies 
that the system of second order differential equations (3.6)-(3.8) is now reduced to first 
order ordinary differential equations which are much easier to deal with. We also point out 
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that the superpotential method is fully equivalent to the standard approach (i.e. solving 
the equations of motion) as long as the solutions for scalar field have monotonic profile, 
as it will be the case in our model with kink-like profiles. It is worth to mention that 
the standard (and straightforward) application of the superpotential method is limited 
to the single scalar-field case since with multi-scalar fields it becomes difficult to handle 
analytically. 

We are interested in the case where the scalar field 4>(y) is given by a kink- like profiles 3 , 

i.e., 



^(y)=5j-^tanh(^(y-y Q )), (3.11) 



a 



where f3 a are the thickness regulators and K a parameterize tensions of the branes in the so 
called brane limit f3 a — > oo (from here on we will consider j3 a = f3, i.e. all the branes have 
equal thickness, although that could be relaxed). As it will be shown in the next section 
the profile (3.11) in the brane limit corresponds to 3-branes with brane-tensions given by 

A Q = (3.12) 

It is important to note that this set-up implies that only positive brane tensions could 
be mimicked by scalar filed configurations, as was also pointed out by DeWolfe et al. [5]. 
Therefore the scalar field can not reproduce the RSI scenario where the IR brane has a 
negative tension. 

In our set-up we consider two kinks corresponding to two thick branes at locations 
y = yi and y = y 2 . They are supposed to mimic two positive-tension branes in the brane 
limit, so the scalar profile (j>(y) could be chosen as follows, 

</>{y) = ^= tanh(/3(y - Vl )) + ^ tanh(/3(y - y 2 )). (3.13) 

We can find the superpotential W((j>) in such a way that it allows a solution of the scalar 
field 4>(y) as in Eq. (3.13). This can be obtained from Eq. (3.10) as, 

dW(cf>) dW(<j>(y)) dy W'(y) 
Hy) = — = -^y—dW) = W (3 } 

W(y) = J {4>\y)fdy + W G , (3.15) 

where Wq is some constant of integration. Deriving the above relation it is assumed that 
(f)(y) is an invertible function of y, so that the superpotential could be written as a function 
of y as 

W(y) =k\ |tanhL% - yi )] - i tanh 3 [/3(y - Vl )\ 

+ k\ |tanh[/3(y - y 2 )\ - ^tanh 3 [/3(2/ - y 2 )]\ + W Q , (3.16) 



3 The profile of the scalar field could be different than the standard kink but the essential concept holds 
for any profile which is monotonic and satisfies equations of motion. 
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where in deriving Eq. (3.16) we assume that the cross term is negligible as far as f3 is large 
and/or the separation "y 2 — yi" between the two thick branes is large such that, 



2/3kik 2 

coslr^(/3(y - yi)) cosh 2 (/3(y - y 2 )) 



d V ,2, a , -TT— TT77T, 7T«0 



After obtaining the superpotential W(y) we find from Eq. (3.10), 

A'(y) = - Ul (tanh(/3(y - yi)) - \ tanh 3 (/?(y - yi)) 

+4 (tanh(/3(y - y 2 )) - I tanh 3 (/3(y - y 2 ))J + W j ■ (3.17) 

The integration constant Wq can be fixed by the requirement that A(y) has a maximum 
at y = uq. The location of maximum with respect to yi i2 will correspond to different 5D 
geometric configurations that we will discuss in Sec. 4. Therefore, we choose the integration 
constant Wo as, 

Wo = - j*i (tanh(/3(y - yi)) - ^ tanh 3 (/3(y - yi)) 

+4 Uanh(/3(y - ifc)) " £ tanh 3 (/3(y - y 2 ))^ | , (3.18) 

such that A'(y ) = 0. 

Now it is straightforward to find the warp factor A(y) by integrating Eq. (3.17) w.r.t. 
y. The result reads, 

+4 { ctfw-ri) - lnCMl,1( * - } + ^p w ° v + Ao ' (3 ' 19) 

where A$ is a constant of integration which can be fixed by the requirement such that 
A(yo) = 0. Note that far away from the thick branes the warp function takes the RS form 
[3, 15] 

A{y ->■ oo) ~ -n\y\, |y| » |yi|, |y 2 |, (3.20) 

where k = 24 ^ /3 (|k| + |k| — Wo). It is easy to see that we get the same behavior, i.e., 
the RS form of A(y), for all values of y in the brane limit when (3 — > oo, which will be 
discussed in detail in the next section. 

It should be mentioned that this set-up reduces to the usual single thick brane [8] if we 
assume that one of the branes is far away from the other brane, say y± = and y 2 — > oo such 
that the second brane can be removed from the set-up and in that case k = qZajs 
Unfortunately in the one-brane case the hierarchy problem remains unsolved, that is why 
we focus on the two-brane scenario as the hierarchy problem is one of the main motivations 
of this work. In our case the second thick brane is located in the AdS space at suitable 
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location 1/2 such that the hierarchy problem could be addressed as in a similar way as it 
was done in the original RSI set-up [3] or in the LR set-up [18]. 

Before closing this section it is instructive to discuss the shape of the scalar potential, 
that is determined by our requirement of having (3.13) together with the ansatz (3.2) as 
solutions of the equations of motion. Having the superpotential W determined, one can, 
using (3.9), find the scalar potential as a function of y. However, since (f)(y) is an invertible 
function, therefore it is also possible to plot the potential V{4>) (3.9), as a function of <j>. 
However, in order to develop some intuition, let us first consider the presence of just one- 
kink profile (2.21) with parameters n and (3. Then the potential V(4>) could be determined 
analytically since one can easily solve equations of motion for A(y), from the invertible 
profile one can find y = y{<f>) and then adopt it, for instance, in (3.6). The results reads: 

Note that the above form of the potential applies, strictly speaking, only for —n/yf]3 < <j) < 
n/^/P since the profile that is required to fulfill equations of motion varies in that range. 
For small field strengths, <f> < M 3//2 , gravitational effects (the second term) are small, so 
that the dominant contribution is just the Mexican-hat potential (the first term). In the 
real case the scalar profile is a sum of two kinks (3.13), therefore its range of variation is 
roughly a "sum" of ranges for two separate kinks. In the absence of gravity each kink is a 
solution of equation of motion for a Mexican-hat type potential, again if a field strength is 
small comparing to the 5D Planck mass 3 / 2 the dominant contribution to the potential is the 
Mexican-hat component. Since kinks profiles vary in between field strengths corresponding 
to the two minima therefore the shapes which we observe in Fig. 1 are roughly "sums" of 
inner parts (only the inner part could be determined) of two Mexican-hat like potentials; 
one centered around —k\/ (2yf]3) and the other one around + K2 / (2\//?) . In reality (with 
gravity) the picture is slightly distorted by the gravity effects that become relevant around 
the external ends of region of variation where 0/A/ 3 / 2 ~ 1. 

In the case of double kink it is not possible to find the potential analytically, so V(4>) 
determined numerically is shown for several choices of the parameters /3 in Fig. 1. Since 
the strength of the profile field varies between —(k\ + K2) /yf]3 and +{k\ + K2) /yf]3 therefore 
the potential V{4>) can also be determined in that region only, which is manifest in Fig. 1. 
Note that in order to trust classical field theory results the scalar field strength (ft must be 
limited by the 5D Planck mass M 3//2 , therefore we conclude that our results are consistent 
if (3 > /c 2 /M 3 . Since in Fig. 1 we assumed ki = K2 = 1 and M = 1 therefore we are limited 
by j3 > 1, so those cases were plotted and then the range of variation of 4> is appropriate. 

4 The brane limit 

In this section we will consider different possible scenarios that could be realized with two 
thick branes and then we will discuss limiting (the brane limit) solutions corresponding to 
thin (singular) branes. To show how scalar field (p is making a brane, we start by looking 
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Figure 1. The potential ^(0) plotted as a function of the scalar field <f> f° r different values of the 
thickness parameter (3 with the same brane tensions Ki = K2 = 1- Hereafter we assume M = 1, 
therefore the field strength is expressed in unites of M 3 / 2 . 



at its action calculated for the profile (3.13), 



dx b ^g~ { -^g MN V M ^N(/> - V{<j>) \ , 



dx 5 ^—g 
+ 



6') 2 + -Lw 2 ' 



6M 3 



6M 3 



cosh (/3(y - yi)) cosh (f3(y - y 2 )) 



k\ tanh(/3(y - yi )) - - tanh 3 (/3(y - Vl )) 



+4 tanh(/3(y - y 2 )) - , tanh 3 (/% - y 2 )) + W Q 



(4.1) 



In the brane limit we have 4 , 



lim i j — — 

/3^oo [ cosh ((3(y - y^) 



-My-Vi) for i=l,2 



such that the scalar action (4.1) can be written as, 



S* = J dx 5 V=g{-^(y - vi) - \i%Kv - 2/2) - A B (y)| , (4.2) 



4 Note that we have corrected the misprint that appeared in [8] in their footnote 3. 
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where A#(y) is a function that generates cosmological constants in various regions of the 
bulk: 



Ae(y) = Jim 

p— >oo 



k\ f tanh(/3(y - yi )) - J tanh 3 (/5(y - yi )) 



6M 3 \ 1 V 3 
-k\ ( tanh(/3(y - j/ 2 )) - ± tanh 3 (/3(y - y 2 ))) + Wo* 



1 f 2 2 1 2 

-Kjsgn(y-yi)) + -Klsgn(y-y 2 )) + W > ■ (4.3) 



Therefore, depending on the choice of the extremum location yo, different values of cos- 
mological constant to the left, in between and to the right of the two branes could be 
generated. 

In what follows we will find analytic solutions for the two positive tension branes in 
the brane limit (/3 — > oo) and numerical results for the corresponding thick brane scenarios. 
In the brane limit, we can write the total action as, 

Sbl = J dx 5 ^ {2M 3 R - M8{ y - yi ) - X 2 5(y - y 2 ) - A B (y)} , (4.4) 

where Ai 2 the respective brane tensions at each brane located at y = y\ and 

y = y 2 and As(y) is the bulk cosmological constant, defined in Eq. (4.3). In the brane 
limit we can obtain the equations of motion from action (4.4) as, 

24M 3 [A') 2 = -A B , (4.5) 
12Af 3 .4" + 24AT 3 (A') 2 = -K B - \ x 8(y - Vl ) - X 2 6(y - y 2 ), (4.6) 

In the brane limit the smooth solution of A'(y) (3.17) will take the following form, 

1 f 2 2 1 

A'(y) = -j^p \ ^\sgn{y - Vl ) + -K 2 2 sgn(y -y 2 ) + W \. (4.7) 

From Eqs. (4.5) and (4.6), one gets (which is also manifested from Eq. (4.7)), 

12M 3 A" = -Ak% - Vl ) - X 2 5(y - y 2 ), (4.8) 

which implies that A" < as A12 > 0, thus alowing for maxima of A(y). Let us denote 
the location of the maxima by yo, that can be chosen anywhere along the extra-dimension, 
such that, 

A'(y ) = 0. (4.9) 
In the brane limit Wo (3.18) takes the form, 

W = - I -n\sgn{yo - yi) + ^lsgn{yo - y 2 ) \ . (4.10) 
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Therefore we obtain (in the brane limit) for the bulk cosmological constant, the following 
result 



1 f 2 2 2 2 

As(y) = -Qj^s \ ^i4sgn(y - y x ) + -nlsgn(y - y 2 ) - -Kjsgn(y - y x ) - -/4sgn(y - y 2 ) 

(4.11) 

It is also important to note that the equation of motion (4.5) implies that the bulk cosmo- 
logical constant is negative leading to anti-de Sitter vacua or in the case where it is zero 
that the corresponding geometry will be Minkowski in that region of space. 

In the following we will consider different cases depending on the location of the ex- 
tremum point yo along the extra-dimension. 

Case-I: We consider the case when the extremum location is on one of the branes, say at 
2/0 = 2/1) we get the analytic results for A'(y) as, 



A'(y) 



24lF A l V < 2/1 

~24lF A i 2/1 < 2/ < 2/2 (4.12) 

~24lP" ( Al + 2A2 ) 2/ > 2/2 



The corresponding bulk cosmological constant A#, in different regions along the extra- 
dimension, reads as, 



> 2 

24M a ' 



2iIF A i 2/1 < 2/ < 2/2 , (4.13) 

2 



"2OT( A 1 + 2A 2) V>V2 

and the values of bulk cosmological constant at the brane locations are A#(yi) = and 
Ab(2/2) = — 12 a/3 (2A2 + A1A2). Note that (4.13) implies in the brane limit correlations 
between the brane tensions Aj and the bulk cosmological constant A#. It is worth to write 
down relations between A#, the slope of the warp function A 1 and the brane tensions Aj in 
the asymptotic regions: 

Ab = f -A! y -> -00 

A' I Ai + 2A 2 y -> +00 v ; 

In the RS2 (with one brane of positive tension A) the corresponding relation is A^/Zc = ±A, 
where k corresponds to ±A'. 

Numerical solutions for the thick branes in the Case I are shown in Fig. 2. This 
configuration is such that the warping function A(y) is positive (negative) to the left (right) 
of the branes, we will see that this scenario will have the normalizable zero- modes of the 
metric tensor perturbations (around the background solution) which correspond to the 4D 
graviton. In the brane limit this set-up is similar to the two positive D3-brane model 
considered by Lykken and Randall in [18]. 

In order to illustrate how are the branes generated, in Fig. 2(a) we show the energy 
density Too = e 2A ((j) 12 /2 + V((f>)) corresponding to the profile. Using equations of motion 
Tqo could be rewritten as e 2A (cf)' 2 — 24M 3 A' 2 ). This form separates two contributions 
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(a) (b) 




Figure 2. Graph (a) shows the shape of the scalar field profile, the corresponding energy density Too 
and also the localized contribution to Too (i-e. e 2A <fr /2 ) for the Case-I with /3 = 3 and Ki = K2 = 1, 
whereas, graph (b) illustrates the shape of the warped factor e 2A , the bulk cosmological constant 
Ab and the potential y[</>(j/)] for the casc-I when the maxima of the warp factor is located on the 
thick brane at y\. 



to Too: local (one that "creates" the branes) e 2A (j)' 2 and non-local oc e 2A A' 2 (one that 
generates the bulk cosmological constants). 

As it is seen from Fig. 2(b) to the left and to the right of the branes the warp factor is 
quickly vanishing, that has been already observed in (3.20). If the branes are sufficiently 
thin (or well separated) then in between them the warping is also nearly exponential so 
that the hierarchy problem could be addressed. We will call the brane located at y\ and 
2/2 as UV and IR branes, respectively. To see the consequences of the warped background 
geometry lets assume that the Higgs field is bounded at the IR brane and its action can 
be written as, 

Sh = ~ J d^x^g [g^d^dvH - m 2 \H\ 2 + A|#| 4 } , (4.15) 

where is the 4D metric induced on the IR brane, g^ = e -2 " 4 ^ 2 - 1 rf" v , with A{jj2) being 
the value of warped factor at the IR brane and m is the 5D Higgs mass parameter (of the 
order of 5D Planck mass). Now the effective 4D action for the Higgs field can be written 
as. 

S H = - J d A x \e 2A ^rf v 3^d v H - m 2 e AA ^\H\ 2 + Ae 4A ^|tf| 4 } , (4.16) 

where we used the fact that, \f—g = e iA ^ y2 \ In order to obtain canonically normalized 
Higgs field, we rescale , H — > e~ A H, such that, 

S H = - J d 4 x \rf v d^d v H - m 2 e 2A{y ^\H\ 2 + A|#| 4 } , 

= - Jd 4 x \rf v d^d v H - n 2 \H\ 2 + A|#| 4 } , (4.17) 

where \i = me A ^ y2 ^ is the effective Higgs mass parameter as viewed on the IR brane. If 
we assume that the fundamental mass scale of the 5D theory is the Planck mass then we 
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can require the value of warped factor at IR brane such that we get the effective 4D Higgs 
mass parameter fi ~ TeV. From Eq. (4.12) we have 

^(y2) = -^3(Ai + A 2 )y 2 , (4.18) 

therefore if 24 ^ f3 (Ai + A 2 ) 2/2 ~ 30 then the hierarchy problem could be solved. Further- 
more, as it will be shown in Sec. 8.1, in this scenario (i.e. yo = yi) there exits a normalizable 
zero-mode which corresponds to the 4D graviton. 

Case-II: Lets now consider the case when the extremum position is in between the two 
thick branes such that yi < yo < U2-, in that case corresponding values of A'{y) (4.7) in the 
brane limit are given by, 



i2iF A i y < vi 

yi<y <V2 (4.19) 

-i2iF A 2 y > yi 



This situation is interesting since the 4D graviton is now normalizable and it is localized 
in between the two positive branes. However in this scenario the hierarchy problem can 
not be solved since there is no warping in between the two branes which is manifest in 
Fig. 3-(a). The two positive branes set-up is now similar to the single brane RS-2 [15]. 
Case-Ill: Now we consider the scenario with extremum located to the left of the left brane 
or to the right of the right brane, so yo < y\ or yo > yi- For the case when the extremum 
lies to the left of yi, in the brane limit, A'(y) (4.7) is given by, 



A'(y) 



y < yi 

~T2iF Al yi <y <y2 , (4.20) 

- T2i^( A i + A 2) y>y2 



In this case we have Minkowski background to left of the brane located at y\ which could 
be called the UV brane, however to the right of this brane we have the warped geometry, 
so that the hierarchy problem could be approached in the same way as it was discussed for 
the case-I. It is worth mentioning that similar geometrical configuration was considered by 
Gregory, Rubakov, Sibiryakov (GRS) [19] with singular branes. The important difference 
is that GRS model have one positive and one negative tension D3-brane while in our case 
the both branes are made out of scalar field which mimic two positive tension branes 
in the brane limit. The numerical results for corresponding thick branes are shown in 
Fig. 3- (h). We would like to comment here that even though this scenario addresses the 
hierarchy problem but it does not have the normalizable 4D graviton, for details see Sec. 
8.1. Similarly the other possibility could be considered when the extremum position is to 
the right of the brane located at y 2 , i.e., with yo > yi, then in the brane limit the A'{y) is 
given by, 



A'(y) 



( Ai + A2 ) y < yt 

12173-A2 yi < y < V2 , (4.21) 

y > y 2 
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(a) (b) 




y 



Figure 3. Graph (a) shows the warped factor e 2A , the bulk cosmological constant A# and the 
potential V[</>(y)] for y± < yg < yi while graph (b) for yo < Vi- Parameters chosen: j3 — 3 and 
«i = k 2 = 1. 

So in this case geometry to the right of the brane at 1/2 is Minkowski. 

In Table 1 we summarize results for the cosmological constant in the brane limit, the 
regions I, II and III are defined as y < yi, y\ < y < y<i and y > y2, respectively. 



Table 1. The bulk cosmological constant Ab when yo 7^ 2/1,2- 



Location of yo 


Region-I 


Region-II 


Region-Ill 


2/1 < yo < 2/2 


1 X 2 





1 X 2 


2/0 < 2/1 







~6m ( Al + A2 ) 


2/0 > 2/2 


-m* ( Al + A2 ) 2 


1 X 2 






In all the above cases, in the brane limit, the A'(y) have discontinuities at the brane 
locations y = y\ and 2/ = 2/2- The discontinuities (or jump) are as follows 

lA ' h = -<m 4 = "12^ " = K (4 ' 23> 

where, [A']i (i = 1,2) is defined as, 

[A'] i = ^m[A'(y i + e)-A'(y i -e)]. 

A jump in A'(y) implies that A"{y) have delta- like singularity, which is consistent with the 
equation of motion (4.6), in fact one can obtain the above jump conditions by integrating 
Eq. (4.6) from yi — e to yi + e and then matching the coefficients of delta functions. 

5 Linearized Einstein equations 

Let us consider fluctuations around the vacuum solution discussed in the previous sections. 
We start by replacing the 5D metric 9MN(x,y) by ffMJv(y) + h M N(x,y), where g~MN{y) is 
the unperturbed background metric, given as 

9fiu = e 2A Vvv, 9^5 = 0, 555 = 1- (5.1) 
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It will be convenient to adopt the Einstein equations in the Ricci form as, 

Rain = Jjj73 T MN, ( 5 - 2 ) 

where 

Tmn = Tmn — -^9mnT a , (5-3) 

where Tmn is the energy-momentum tensor given in Eq. (3.5) while = — |(V(/>) 2 — 5V(4>) 
is its trace. So we get, 

2 

T M n = Vm0Vat0 + ^9mnV((P). (5.4) 

The perturbations in the Tmn will correspond to fluctuations of the scalar field 4>(y) as 
4>(x,y) = 4>{y) + (p(x,y) and of the metric gMN(x,y) = 9m N (y) + h M N(x,y). These 
perturbations can be calculated order by order in perturbation expansion as, 

Tmn = T$ N + t mn ~^ > 

where ellipses correspond to the higher order fluctuations in tp{x,y) and hMN(x,y)- The 
zeroth and the first order terms are as follows, 

f^ = 2 -e^V^), 4°W 2 + ^), 2*?=0. (5.5) 

^ = ^(e 2 ^^^ + ^)v)> (5-6) 

f « = 20V' + \*^-<P + j^W). f « = + \h^V{ct>). (5.7) 

Using results from the Appendix A one can find the following explicit expressions for 
components of the Ricci tensor in the zeroth and first order, 

= - e 2A (A" + AA' 2 ) r^, 4? = "4 (A" + A' 2 ) , fljg = 0, (5.8) 
R$ = ~ \d^d u h 55 + e 2A Vflv {A" + AA> 2 ) h 55 + ±e 2A V ^A'ti 55 + \ [d^ + 

+ A' + £U M s) - \e~ 2A Uh^ v + V 2 V CT (a M 5 p /w + dvdphp, - d^8 v h pa ) 

- l -h% - ±A' VtiU r)0°ti p<T - A' 2 (2^ - + A'v»T9^, (5-9) 

<5 =^ e ~ 2A ^ - d^) - e-^A'rT (5 P V " ^ V) 

+ ^'^/iss - ^e" 2A {Uh^ - rf°WM - {A" + 4A' 2 ) h^, (5.10) 

(5.11) 

where □ is the 4D d'Alambertian operator, i.e., □ = tfdJd v . 
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Having all the components of the Ricci tensor (5.9)-(5.11) and Tmn (5.6)-(5.7), one 
can write down the equations of motion for the metric fluctuations liMN(x,y): 

(H : - \d p 8 u h 55 + e 2 \ pv {A" + 4A' 2 ) h 55 + ^e 2A Vflu A'ti 55 + \ (8^5 + ^5) 

+ a' (d„K 5 + d u h^) - \e~ 2A uh pv + V 2 V" (d^dphn + Qudph^ - dpd v h pa ) 



\K V - ^A'rj^h'^ - A' 2 (2/v - ri^rfhpa) + A' V ^d p h a5 



1 2/ "e 2 V^y + ^)V)> (5-12) 



4M 3 3 



0*5) : Je" 2 (3,/^ - dph' pa ) - e- 2A A'fT (d p h pa - d p h pa ) + U%h 5 



(55) 



- ^e~ 2j4 (D^ 5 - rT-V^) = (5.13) 
e -2A ^y*^ + A"r) pCT h pa - \r} pa h" p )j - \e~ 2A Uh^ + 2^5 + e" 2 V^Xs 

2^V + |^^ + ^5 5 n0)V (5-14) 



4M 3 V 3 <9</> 3 



Additionally, we also have the equation of motion of the scalar field eft (3.4) in the first 
order in the fluctuations hMN(x,y) and ip(x,y) as, 



e~ 2A Oip + tff + 4AV - + ^' (e" 2 ^)' - ^'hf - (<j>" + 4A'^) h\ = 0, (5.15) 



— = — id 4- 

2^ v ; 2 

where ft = rj^h^. 



6 SVT decomposition of perturbations and gauge choice 

In all expressions derived in the previous section the perturbation Km n represents the most 
general symmetric small fluctuation around the background. Due to the symmetries (4D 
Poincare invariance) of the background metric and the energy-momentum tensor, we can 
decompose the perturbations Kmn into scalars, vectors and tensors as follows 

V = e2A ~ M^E + d^G v + 3 V G P + H/jv] , (6.1) 

hps = dpB + (6.2) 
h 55 = 2 X , (6.3) 

where ip, Xi B an d E are scalars, whereas, C M and G p are diver genceless vectors and 

is the transverse and traceless tensor, for details see the Appendix B. The perturbation 

modes are functions of x and y coordinates. 

Now we will consider the coordinate/gauge transformations and then we will turn to 
the question of choosing the appropriate gauge in order to eliminate artifacts of the freedom 
of choosing a reference frame. Lets consider the following coordinate transformation, 

x M = x M -^ M , (6.4) 
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where the £ Ajf is an infinitesimally small function of space time, i.e., \£, M \ << \x M \ and 
£, M = (C^ 1 ) £ 5 ) with being a 4D vector and £ 5 a scalar change in the 5th coordinate y. In 
order to write down corresponding gauge transformations of the decomposed scalars, vector 
and tensor modes, it is useful to decompose also the 4D vector £ M into the divergenceless 
vector £j_ and gradient of the scalar £11, i.e., 

e = e ± +o^ a^ = o. (6.5) 

It is easy to show that the change in metric perturbation Hmn corresponding to (6.4) reads 

hMN = h M N + 8h M N, (6.6) 

with 

ShMN = Vm£n + Vjv^m, (6.7) 

where as usual Vm is the 5D covariant derivative, see the Appendix A. The explicit form 
of the components of SHmn are given by, 

*V = d ^ + + 2d ^\\ + 2A'e 2A r]^ 5 , (6.8) 

Sh,* = dpts + + " 2A 'Z±» ~ 2A%^ , (6.9) 
5h 55 = 2& (6.10) 

The above transformations of the metric perturbation Hmn induce the corresponding trans- 
formations of the metric perturbation components defined by Eqs. (6.1)-(6.3) as, 

?P = ^-A'f, E = E-e~ 2A i h (6.11) 

X = X + C' 5 , B = B + i[+e-2A'i h (6.12) 

C^ = C^ + t'^-2A't^ G^ = G^ + e~ 2A ^ (6.13) 

whereas, is unaffected by the coordinate transformations. Similarly, the gauge trans- 
formation of the scalar field perturbation <p can be easily obtained as, 

(p = tp + 5(p = (p + (f>'f. (6.14) 

Similarly, the gauge transformations of the energy momentum tensor can be written as, 



where, 



^MN ~ ^MN + ^MNi (6.15) 



6fti> N = fj°\v N Z A + tSv m ? b + V C f$> N f- (6-16) 



For completeness we show also the explicit form of the transformation for the {^v) com- 
ponents in terms of decomposed £: 

Sf<$ = \e 2A V{4>) {d^ + d^±, + 28^ + 2A'e 2 %M) + j^V^V^ 

(6.17) 



- 20 - 



where we have used the zeroth order value of T^J (5.5). Similarly, one can derive trans- 
formations for the other components but this will serve our purpose. 

Now we turn our attention towards the issue of choosing a gauge. It proves to be con- 
venient to adopt the so-called longitudinal or Newtonian gauge defined by the conditions: 
B = E = in the scalar and = in the vector sector. It is important to note that, 
indeed, one can always choose the gauge parameters such that the gauge conditions are 
satisfied, i.e., 

^(x,y) = e 2A E, (6.18) 

?(x,v) = -B-f A {e- 2A E)', (6.19) 
Sxn(x,y) = -e 2A G„. (6.20) 

Here one can note that this choice of gauge fixing completely fixes the gauge and so that 
there is no residual gauge freedom left. 



7 Equations of motion for the SVT modes 

In this section we will derive equations of motion for perturbations of the metric and 
the scalar field. We are going to adopt the decomposition (6.1)-(6.3) for Hmn- As it is 
shown in the Appendix B (B.15)-(B.21) the SVT modes decouple, therefore in the following 
subsections we discuss them separately case-by-case. 



7.1 Scalar perturbations 

Scalar perturbations contribute to the metric as follows 

ds 2 = e 2A [(1 - 2tp) ft,,,, - 2d li d v E\ dx^dx u + d^Bdx^dy + (1 + 2 X ) dy 2 , (7.1) 

The scalar modes appearing here are not gauge invariant, i.e., their values are affected by 
the choice of different coordinates. It is therefore instructive to either work with the gauge 
invariant quantities or choose a suitable gauge such that the ambiguities related to the 
coordinate transformations can be removed. Here we choose the longitudinal gauge such 
that the gauge freedom is fixed completely, as discussed earlier in Sec. 6. Therefore, for 
the scalar modes of perturbation, we have B = E = 5 . 

In the longitudinal gauge the perturbed metric (7.1) is of the form, 

ds 2 = e 2A (1 - 2i/>) rj^dx^dx" + (1 + 2 X ) dy 2 . (7.2) 



5 We suppress the signs hereafter, as it is clear that we are referring the the modes in the new reference 
frame. 
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Adopting the general results from the Appendix B we find the following form of the lin- 
earized field equations for the scalar modes, 



(/XI/) 

(/x5) 
(55) 



2 (A" + 4A' 2 ) X + A' X ' + e- 2A DV + 8A'ip' + 

1 



+ d„d v (2^ - x) 



6M 3 



2A dV{4>) 



1 



4 + 2A>') + 4A' X ' - e~ 2A Ux 



(7.3) 
(7.4) 



4M 3 



(7.5) 



One can notice from Eq. (7.3) that the absence of the d^d v term on the right hand side 
implies d^d u (2i/j — x) = so that x = 2^ + c (y)- Where c(y) is a y-dependent constant of 
integration which can be fixed by the requirement that at 4D infinities Xi ift — > 0, therefore 
c(y) = 0. 

When x = 2^> the equation of motion for the scalar field fluctuation (5.15) simplifies, 

3 2 V(0) 



-<p - 6</>V - 4 (</>" + ^ = 0. 



(7.6) 



It is important to note that, as usually in such cases, the equations of motions (7.3)-(7.6) 
are not independent. Adopting the relation x = 2V> and the background equations of 
motion one derive the following equation that we will use instead of (7.3) and (7.5) 

1 



3^" + 6^>' - 3e- 2A \DiP 



2M 3 



6y. 



(7.7) 



Hence, Eqs. (7.4), (7.6) and (7.7) complete the set of linearized equations for the scalar 
modes. In the subsequent section we will use these equations to study stability of scalar 
field perturbations. 



7.2 Vector perturbations 

We can write down the metric for the vector perturbations as, 



e 2A (vuv + d^G v + d u G^) dx»dx u + C^dx^dy + dy\ 



(7.8) 



where and are divergenceless vectors defined in Eqs. (6.1) and (6.2). Adopting 
the general results from the Appendix B we find the following form of the linearized field 
equations for the vector modes 



0„ 



C' v + 2A'C U + -Ae 2A A'G' v - e 2A G" 



-2A 



UC V - DG'„ 



0. 
0. 



(7.9) 
(7.10) 



Since we are working in the gauge where G^ = so the equations of motion for the vector 
modes of the metric perturbations read 



UC V = 0, (C' u + 2A'C V ) = 0. 



(7.11) 
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7.3 Tensor perturbations 

The tensor metric perturbation (2.2) can be written as, 

ds 2 = e 2A M( V ^ + H^ u )dx»dx u + dy 2 , (7.12) 

where, H^ u = H fJjL/ (x,y) is the tensor fluctuation as defined in (6.1). Adopting the general 
results from the Appendix B we find the following form of the linearized field equations for 
the tensor modes 

(d! + 4A'd 5 + e~ 2A U) = 0. (7.13) 

The zero-mode solution (corresponding to = 0) of the above equation should rep- 

resent the 4D graviton while the non-zero modes are the Kaluza-Klein (KK) graviton 
excitations. 

The issue of the stability of the background solution (in terms of the SVT components) 
will be discussed in the next section. 



8 Stability of the solutions 
8.1 Tensor perturbations 

In order to gain more intuition and understanding of the graviton equation of motion (7.13), 
it is convenient to change the variables such that we get rid of the exponential factor in 
front of the d'Alambertian and the single derivative term with A', so that we convert the 
above equation into the standard Schrodinger like form. We can achieve this in two steps; 
first by changing coordinates such that the metric becomes conformally flat: 

ds 2 = e 2A ^ {g^dx»dx u + dz 2 ) , (8.1) 

where dy = e A dz. In the new coordinates the Eq. (7.13) takes the form 

(d 2 + 3A(z)8 z + □) = 0, (8.2) 

where dot over A represent the derivative with respect to z coordinate. Now we can 
perform the second step removing the single derivative term in (8.2) by making the following 
redefinition of the graviton field 

H^(x, z) = e- 3A / 2 H^ u (x, z). (8.3) 

So the Eq. (8.2) will take the form of the Schrodinger equation, 

(d 2 z - \\ 2 {z) - \A(z) + D) H, u (x, z) = 0. (8.4) 

We can further split the H^ u (x,z) into i^,(x, z) = H fiV (x)H(z), where H fJiV (x) = e tpx is a 
^-independent plane wave such that nH flu (x) = m 2 H^ v (x), with — p 2 = m 2 being the 4D 
KK mass of the fluctuation. Then the above equation takes the form, 

[d 2 z - \A 2 {z) - \A{z) + m 2 ) H{z) = 0, (8.5) 

[-d 2 + U{z)] H{z) = m 2 H{z), (8.6) 
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where U(z) is the potential, 



U(z) = \A\z) + (8.7) 



Note that we can write the Schrodinger-like equation (8.6) in super symmetric quantum 
mechanics form as, 

Q^QH = f-d z - lA\ (d z ~\Ah = m 2 H. (8.8) 

The zero mode (m 2 = 0) profile, Hq(z), corresponds to the graviton in the 4D effective 
theory. The stability with respect to the tensor fluctuations of the background solution 
is guaranteed by the positivity of the operator Q)Q in the supersymmetric version of the 
equation of motion (8.8) as it forbids the existence of tachyonic modes with negative mass 2 , 
m 2 < 6 . So, in that case, the perturbation is not growing in time. Hq(z) could be found 
by noticing that the annihilation operator Q should vanish acting on Hq, 

QH = (d z - H = 0, (8.9) 



2 

which implies that, 

H (z) = ei A{z \ (8.10) 



For massive KK modes we need to solve the Eq. (8.6) with m 2 ^ 0. For large z the potential 
U (z) goes to zero for the case (i) and (ii) as shown in Fig. 4(a) and (6), so Eq. (8.6) reduces 
to one dimensional Klein-Gordon (KG) equation, i.e., 

(8 2 + m 2 ) H m {z) = 0. (8.11) 

Therefore in the large z limit, we expect, 

H m (z) « c\ cos(mz) + C2 s'm(mz), (8-12) 

where c\ and C2 are constants. Therefore the massive KK modes are plane wave normaliz- 
able and we have a continuum spectrum of KK states for the cases {%) and (ii) discussed in 
Sec. 4. Before closing this subsection we will briefly comment on the effective 4D gravity. 
We are going to estimate the effective 4D Plank mass and discuss the localization of the 
zero-mode of the perturbation and then corrections to the Newton's potential due to the 
massive KK modes. To calculate the 4D Plank mass it is important to note that Eq. (8.4) 
only involves 2nd derivatives of the metric perturbation H fJjU (x,z) which is related to the 
fact that in the action these fluctuations have the following canonical kinetic term 

S » M 3 J d A xdzd M H^(x, z)d M H tMV (x, z) + ■ ■ ■ , (8.13) 



6 Since J dz(QH) 2 + H QH ft" = m 2 f dzH 2 , therefore in order to guarantee m 2 > the boundary term 
must vanish or be positive. Note that for the zero mode QHo = 0, so indeed the boundary term disappears. 
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(a) (b) 




Figure 4. These graphs illustrate the shape of the quantum mechanics potential U(z) in solid 
(black) for all the for scenarios that we have considered in Sec. 4 and the corresponding shape of 
the zero-mode (4D graviton) in dashed (red) curve. Parameters chosen: /3 = 5, Ki = 3 and K2 = 1. 



where the indices are contracted with the 5D Minkowski metric T)mn- Since we have 
H^ u (x, z) = H flv (x)H(z) therefore (8.13) could be rewritten as follows 

S » M 3 J dzH 2 (z) J d A xd a H flu (x)d a H f * u (x) + ■■■ , (8.14) 

from this we can read out the effective 4D linearized gravity as, 

S « M 2 Pl j d A xd a H l , u (x)d a H flu (x) + • • • , (8.15) 

where Mpi is the effective 4D Planck mass, i.e., 

M Pl = M 3 f dzH 2 (z), (8.16) 



where H{z) satisfies the supersymmetric quantum mechanic equation (8.8) for m 2 = 0. In 
order to reproduce the standard 4D General Relativity, M pl must be finite, in other words 
Hq(z) must be normalizable. It is easy to see from (8.10) that indeed Hq(z) is normalizable 
for the cases (i) and (ii) for which the warp function A(y) posses the following asymptotic 
behavior (see Eqs. (4.12) and (4.21)), 

A'(y) < as y -> 00, (8.17) 
A'(y) > as y -»• -00.. (8.18) 
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The above implies that 

J dzHl{z) = j dze 3A ^ = J dye 2A ^ < oo, (8.19) 

therefore Ho(z) is normalizable (see also Fig. 4 (a) and (b)) and Mp t is finite for the cases 
(i) and (ii). The situation for the case (to) is far more complicated, as there neither we 
have the finite 4D effect Planck mass nor we have a normalizable zero-mode (see Fig. 4 (c) 
and (d)). However, as it was pointed out for similar singular brane set-up (GRS [19]), the 
effective 4D gravity on the brane can be reproduced and we could have the quasi-localized 
gravity [23, 37-39]. We are not going to discuss that case any farther. 

In order to build some intuition concerning possible size of corrections to the Newton's 
law due to an exchange of massive KK modes, we will briefly consider a simple case with 
two point like sources mi and ni2 located on the 4D slice at z = z\ = and at z = Z2 
in the 5th dimension, so at locations of perspective branes in the brane limit. With these 
assumptions, the corrections to the Newton's law for the case (i) (the thick brane version of 
Lykken-Randall model) could be easily obtained. At the UV brane (z\ = 0) the potential 
(after integrating over graviton KK modes) is modified as follows 



r 



-mr 



U uv (z) « G N ^ + ±. r^^— iUO) 2 , 



where Gn = (87rMpj) _1 is the 4D Newton's constant, H m (0) is the value of the graviton 
wave function at the UV brane and r is the distance between the point sources on the 4D 
slice located at z = z± = 0. Now, if sources are located at z = Z2 the gravitational potential 
is modified as follows 

Ui R (z) G n h -r-7o / dm e OA(Z2l H m (z 2 ) , (8.21) 

r M- 3 J r 

where H m {z2) is the value of the graviton wave function at the IR brane. For singular 
branes the exact analysis performed in various set-ups show [15, 18, 19, 23, 26, 37-39] that 
corrections to the Newton's law are small and the usual 4D gravity is restored on IR branes, 
therefore one could expect that similar conclusions hold also for our thick brane scenarios. 
In fact, as it was pointed out by Csaki et al. [7] also in thick brane scenario corrections due 
to massive KK modes do not induce any harmful effects for the 4D effective gravity. The 
detailed and rigorous study of these issues for our two thick branes, is beyond the scope of 
this paper and will be considered somewhere else. 

8.2 Scalar perturbations 

The linearized field equations corresponding to the scalar modes of the perturbation are 
given by Eqs. (7.4), (7.6) and (7.7), i.e., 

GA'drf + 38^' = jj^t'd^, (8.22) 



i 



e-^Op + / + 4AV - ~ 60V - 4^^V = 0, (8.23) 



^" + 2AW - e- 2A DV = g^V- (8-24) 
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One can integrate Eq. (8.22) over ^-coordinates and get the following equation, 

1 



6AV + 3-0' 



4M 3 



.25) 



where we have put the y-dependent integration constant to zero by the requirement that 
the perturbations vanish at 4D infinities. It is more convenient to use the conformal frame 
where the metric can be written as in Eq. (8.1) such that dy = e A ^dz. Hence, in the new 
coordinates our equations of motion (8.23)-(8.25) take the following form, 



□93 + (p + 3Aip — e 



2 A d 2 V{ct>) 



^-6^-4e-^ = 0, 



4> + Aip — Dip 
2Aip + ip 



1 



6M 3r 
1 

12M 3 "' 



First we solve Eq. (8.28) with respect to ip and calculate (p as, 

12M 3 u 



r 



2A\p + 2Aip + ip)cp- ( 2Aip + ip)<p 



and then use it in (8.27), so that we obtain an equation only for ip, 



(8.26) 
(8.27) 
(8.28) 

(8.29) 



ip + [3A - 2^ V + 4A - 4A-^ + □ W) = 0. 



(8.30) 



To convert this equation into the Schrodinger form it is instructive to remove the first 
derivative terms of the perturbation tp, to do so we redefine the scalar perturbations as, 



ip(x,z)=e ^ A ^<j)ip(x,z). 



(8.31) 



Then the linearized field equation for the scalar perturbation ip takes the following form, 

/ .. x 2 



-ip + 



lA 2 --A + At + 2 



ip = Dip. 



(8.32) 



We can further decompose the ip(x,z) into tp(x,z) = ip(x)ip(z), where ip(x) = e ipx is a 
^-independent plane wave such that O0(a;) = m 2 ip(x), with — p 2 = m 2 being the 4D KK 
mass of the fluctuation. So, with this field decomposition Eq. (8.32) can be written as, 



-4>(z) + 



»A 2 -°-A + A% + 2 
4 2^6 



ip(z) = m 2 ip(z). 



5.33) 



The properties of this equation has also been explored in the past in the context of stability 
and dynamics of radion in [9, 24, 25] and [26]. To develop some intuition concerning this 
equations it is convenient to rewrite it in supersymmetric quantum mechanics form. For 
this purpose we introduce an auxiliary function a(z) defined by 



a(z) 



e§ A &d>(z) 



A(z) 



(8.34) 
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Now we can write the potential of the above equation in the following form, 



-A 



-A + A-V + 2 ~ 



a{z)d\ 



a(z) 



uj 2 (z)-6j(z), (8.35) 



where u){z) = 
mechanics form as, 



Then we can rewrite the Eq. (8.33) in a super symmetric quantum 



-d 2 tp + (oj 2 {z) — io{z)) tp = m 2 -0 



m 2 ip, 



where the operator A^ and A are defined as, 
A* = (-d(z) + u(z)) , 



A=(d(z)+u(z)). 



.36) 



.37) 



The above supersymmetric form of the scalar perturbation equation (8.36) guarantee that 
there is no solution for with m 2 < 0, hence the fluctuation -0 can not destabilize the 
background solution. The zero-mode for the scalar perturbation ijj(z) can be obtained from 
(8.36) as, 

i A(z) 



a(z) 



(8.38) 



which is plotted as a function of z in Fig. 5. It is important to note that the zero-mode 
for the scalar perturbation ip{z) is not normalizable. As one can easily see from the above 



expression, ipQ diverges when z — > oo (then 



and also es^ 2 ' — > for the case (i) 



and (ii)), this behavior is illustrated in Fig. 5. In fact, also tpo(x, z) 



e 2 



^(x)e~ 3A ^ A{z) is not normalizable. 

Although our main concern here is to verify the stability, nevertheless it is worth to 
check the behavior of the potential at z — > =Loo in order to see whether there is a mass gap 
in the spectrum of scalar modes. From (8.35) one can easily find the explicit form of the 
potential as a function of z 



U^z) 



1a" + -A' 2 + 2A'^- + 2[^- ) 



y=y{z) 



where y as a function of z could be determined from 

ry 



- A ^dy' = z(y)-z(y ). 



.39) 



5.40) 



If we limit ourself to the large y region and the integration constant yo is large enough, we 
can use the asymptotic behavior of A = A(y) as in (3.20) then we find 

1 



~ — ln(KZ + const.) 



where k = kjttz (f + ~ Wo)- From (8.39) we find that \\m z ^ ¥ ± 00 U 1 p{z) 
therefore we conclude that the spectrum is continuous starting at m 2 = 0. 



3.41) 
\ 



7 Of course, X\vci v ^,± 00 U^[z(y)] = 0, as well. 
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Figure 5. These graphs illustrate the shape of the quantum mechanics potential U^(z) (8.35) in 
solid (black) for all the for scenarios that we have considered in Sec. 4 and the corresponding shape 
of the scalar zero-mode in dashed (red) curve. The parameters are same as in Fig. 4. 



It is worth to comment on another possible zero mode solution. The theory that we 
are discussing here is invariant with respect to a shift along the extra dimension: y — > y + e, 
therefore if a given metric gMN(x,y) and a scalar field cj>(x,y) are solutions of equations 
of motion, then so are gMN(x,y + e) and <f>(x,y + e). Expanding them around e = one 
obtains 

9MN(x,y + e) = 9MN{x,y) + gMN{x,y)'e-\ (8.42) 

4>(x, V + e) = 4>(x, y) + (f)'(x, y)e H , 

where ellipsis stand for higher powers in e. Since gMN(x,y + e) and (j)(x,y + e) and also 
9mn(x, y) and cj)(x, y) satisfy the equations of motion, therefore qmn(x, y)' and 4>'(y) satisfy 
linearized equations of motion. In our parameterizations of the perturbations, (6.1)-(6.3), 
that corresponds to 

il>(x,y) = -A'(y), <p(x,y) = <f>'(y) and B = E = X = (8.43) 

As ijj and ip given by (8.43) correspond to modifications of the field configuration (that 
satisfies the equations of motion) along the symmetry directions therefore it is supposed to 
be a zero mode. Indeed, it could be verified explicitly that ip and (p given by (8.43) satisfy 
linearized Einstein equations (7.3)-(7.5) together with the scalar field equation of motion 
(5.15). It should be emphasized that in this case the relation d^dy (2ijj — \) = does not 
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hold by the virtue of 2ip — x = 0, but by the fact that tp(x, y) is ^-independent while x = 0- 
We will not consider those modes any more since they do not depend on x and therefore 
can not be localized in 4D. 

8.3 Vector perturbations 

The field equation obtained for the transverse vector mode of the perturbation, after inte- 
grating Eq. (7.11) w.r.t. x-coordinate, is, 

DCp = 0, C* + 2A'C^ = 0, (8.44) 

where we have set the integration constant to zero by using the fact that perturbations 
should be localized in 4D so that they do vanish far away from sources. It is more intuitive 
to write the vector perturbation in the conformal coordinates so that the results can be 
interpreted easily. Therefore, in the conformal frame the equations of motion for the vector 
modes of the perturbation take the form, 

UCp = 0, + zAc^ = 0. (8.45) 

One can immediately notice from Eqs. (8.45) that the vector modes of perturbations are 
massless. 

Since the Eq. (8.45) is first order in z-derivatives so it can not be put into an elegant 
Schrodiger like form as for the case of tensor and scalar modes. Therefore to see if these 
modes are localized or not we have to find canonical normal modes of these perturbations 
from the second order perturbation of the action [24], the result reads: 

5 2 S V = J d 5 x^ (r]^d^C a d v C a ) , (8.46) 

where, Cn = e^ A C fl corresponds to the canonical normal mode. From Eq. (8.45), one finds 
that Cu(x,z) = C(x)e~ 3A . So the canonical normal zero-mode of the vector perturbation 
can be given as, 

C M = e-^C^x), (8.47) 

where satisfies the equation OC^(x) = 0. Recall that from the requirement of repro- 
ducing the General Relativity at low energies we had 

M| = M 3 J dze 3Aiz \ (8.48) 

Therefore the canonical normal vector modes can not be localized since the integral J dze~ 3A ^ 
must be divergent (as a consequence of the finiteness of the 4D Planck mass). Hence, the 
vector modes of the perturbation are not localized and therefore they do not affect issue 
of stability. 
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9 Conclusions 



Five dimensional Randall-Sundrum like models offer an elegant and simple solution to the 
hierarchy problem. The standard formulation of those models assumes the presence of 
infinitesimally thin branes embedded in five dimensional space time. Usually, one of the 
branes has a negative tension. This work was motivated by a desire to avoid infinitesimally 
thin branes and instead to model them by physical objects, e.g. background profiles of a 
scalar field. 

We have shown that even in the presence of non-minimal scalar-gravity coupling it is 
not possible to mimic a negative tension brane. In an attempt to construct a model that is 
periodic in the extra dimension we have derived a generalization of the Gibbons-Kallosh- 
Linde sum rule that holds also if the scalar field couples non-canonically to the Ricci scalar. 
It turned out that even in that case periodicity forbids any non-trivial scalar field profile 
along the extra dimension. Therefore we have focused on non-compact extra dimensions. 
In order to have a chance to address the hierarchy problem, the scalar background that we 
introduce is composed of two kink-like profiles. This set-up in the brane limit corresponds 
to a model with two thin branes both having a positive tension. Various possible cases, 
depending on the location of the maximum of the warp function has been considered; the 
most attractive option turned out to be the one with the maximum located on the top of 
one of the thick branes. 

Stability of the background solution was discussed in details and was verified in the 
presence of the most general perturbations of the metric and the scalar field. 
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A Conventions 

In this paper we used the metric signature as — h + + + and in our conventions the capital 
roman indices represent 5D objects, i.e., M, N, ■ ■ ■ = 0, 1, 2, 3, 5, whereas, the Greek indices 
label four-dimensional (4D) objects, i.e., fx, v, • • ■ = 0, 1, 2, 3. Some of the most frequently 
used quantities are summarized in this appendix. In our conventions the definition of 
5D covariant derivatives, acting on contravariant and covariant vectors are, respectively, 
V M V N = d M V N +T^ A V A and V M V N = d M V N -r^ N V A - The definition of the covariant 
derivative of a second rank tensor is, 

V A T*f = d A T»f + T%T» - T C AN T^. (A.l) 
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The 5D d'Alambertian operator V is defined as, 



V 2 = V M V M = -^d M ^gg MN d N . (A.2) 



-g 

The metric is perturbed as, gMN = 9m n + h-MN, whereas, the inverse of metric per- 
turbation is h MN = —g MA fiABg BN ■ The unperturbed (background) metric satisfies the 
following ansatz: 

ds 2 = e 2A{y) Vflu dx"dx u + dy 2 , (A.3) 

where the warp function A(y) is only a function of the extra-spatial coordinate y. 

The Christoffel symbols or the affine connections take the following form in terms of 
the linear perturbation of the metric, 

^mn = 2^ AB \P M ^ NB + ^n^mb — dshMN — ZhBC'F mn] i (A-4) 

where the barred quantities, i.e., g and f, are unperturbed. The only non-vanishing com- 
ponents of the unperturbed Christoffel's symbol are, 

f£ 5 = f t = A'SS, f% = -A>e 2A rj, u . (A.5) 

The non-zero components of the perturbed Christoffel's symbols to the first order in the 
perturbation are, 



T p = 




(A.6) 


T-l/i _ 

1 1/5 — 




(A.7) 


r 5 = 

fJ,U 


2 [d/^5 + c^fytS - h'^ + 2A'e 2j4 /i 55 ^^] , 


(A.8) 


T-l/i _ 

1 55 _ 


ie" 2 V" [2^5 - 9^55] , 


(A.9) 


r 5 - 


- [^/iss - 2A'v] , 


(A.10) 


r 5 - 

1 55 — 




(A.11) 



The following relation for Christoffel's symbols proved to be very useful: 

r%v = \ds [e~ 2A K + hi) • (A.12) 

The explicit form of 5D d'Alambertian operator acting upon an unperturbed (4>(y)) and a 
perturbed quantity, say (p(x,y) = 4>(y) + (p(x,y) are, 

V 2 ^y)=4>"(y) + 4A f (y)fi(y), (A.13) 
V 2 0(x, y) =e~ 2A U^ + / + 4AV + ^' (e' 2A r,^ V)' " ^ W " + 4A'0') ^ 

-^V"(«9 M ^ 5 + ^V), (A.14) 

where □ = r} pu d^d u is the 4D d'Alambertian operator. The Ricci tensor in the first order 
in the perturbation can be written as, 

RmN = ®A^MN ~ 9m^AN + ^AB^MN + f AB^MN ~ ^MB^AN ~ ^MB^AN- (A. 15) 
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B Decomposition of the metric perturbation 



We decompose the background perturbations Hmn hito scalars, vectors and tensors as 
((6.1)-(6.3)), 

V = e 2A [-*hv» ~ ^d v E + d„G v + d y G^ + H„ v ] , (B.l) 

h flS =d li B + C fl , (B.2) 

h 55 = 2 X , (B.3) 

where tp, x, B and E are scalars, whereas, and are diver genceless vectors and 
is the transverse and traceless tensor, i.e., 



d"Cp = d"G„ = 0, d^H^ = H% = 0. (B.4) 

Let us discuss the uniqueness of the above decomposition. It is easy to see that B is 
determined by as follows 

UB = d*h&. (B.5) 

Therefore shifting B by a solution the homogeneous equation DA = leads to another 
allowed solution of (B.5) 8 . In order to specify the solution of DA = one has to fix initial 
conditions, that can be done e.g. by specifying X(t, x, y) and dt\(t, x, y) at a given time. 
Hereafter we are going to assume that at a certain time t = to that is far enough in the 
past both X(t, x, y) = and dt\(t, x, y) = 0. That assumption is physically well motivated 
as there is no reason to observe any perturbations at the very beginning and implies that 
the only solution of DA = is in fact A = 0. Therefore the decomposition (B.2) is unique. 
Similar strategy could be adopted to show uniqueness of the decomposition of h^ v provided 
appropriate initial conditions are adopted. We start by determining £ as a solution of the 
following equation that is implied by (B.l): 

r? E = \e~ 2A (hjh» - v) • (b.6) 

Having E determined (with appropriate initial conditions that ensures uniqueness) one can 
find tp solving 

V> = --e- 2A h» - -DE. (B.7) 

Then G^ is a solution of 

□G> = e~ 2A d u V + 2^ (V + UE). (B.8) 

8 Another way of seeing the same freedom in determining B and C M is to notice that a shift B — > B + A 
can be compensated by an appropriate change of C M , C M — ¥ C M — <9 M A. Requiring DA = 0, guaranties that 
C M remains divergenceless. 
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Now we can write down the first order Einstein equations in terms of the scalar, vector 
and tensor (SVT) components defined in (B.1)-(B.3) as, 



(fits) 



0*5) 
(55) 



ip + 8A'ip' + 



2 (A" + AA' 2 ) x + A' x ' + 2 (a" + AA' 2 + \e~ 2A D^ 

+ 8^d u B' + 2A'd^8 u B + A' V ^DB + i (fyC£ + d^) + A' (d„C v + 
+ e 2A [2 (A" + 4A' 2 + □) d^d v E + n^A'UEf + lA'd^Ef + fy&E"] 

+ a M ft, W -x)- [a" + 4A' 2 + ±nj H, u - 2e 2A A'H' tll/ - \e 2A H% 



1 2 



,2-4 



4M 3 3 



V^^jrr-V + ^(0) (-2^ - 2c>^£ + ^G, + d v G» + fl^) 



4 + 2AV) + 4AV - e" 2A Dx + □ + 2A'£') + e~ 2A UB' 



4M 3 



(B.9) 
(B.10) 

(B.ll) 



Adopting the background equations of motion the above equations could be simplified as 
follows: 



(jus) 



( M 5) 
(55) 



e 2A V^ 



1 



2 (A" + 4A' 2 ) x + A'x' + e~ 2A D^ + BA'if/ + V" + e^A'DS + A'UE' 
2?p - X + B' + 2A'.B + e 2A (2QE + 4A'£' + E") 
C' + 2A'C V - 4e 2A A'G'„ - e 2A G',' 



C'„ + 2A'C U - Ae 2A A'G'„ - e 2A G'' 



1 (DH, U + 4e 2A A'H'^ + e 2A H%) 



1 2 



2.4 



4M 3 3 



a^(</>) 



4 (-0" + 2AV) + 4AV - e~ 2A Ux + □ + 2A' E') + e^DS 

1 



4M 3 



(B.12) 
(B.13) 

(B.14) 



Now comparing the coefficients of r]^ u , d^d u , d v and the tensors on both sides we get from 
(/uz/) components the following equations of motion for the scalar, vector and tensor modes 
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of the perturbations 



2 (A" + AA' 2 ) x + A' X ' + e~ 2A DiP + 8A'^' + < + e - 2A ATJ5 + A'Dtf = J_-^M, 

4M 3 3 d$ 

(B.15) 
(B.16) 



2^ - X + B' + 2A'B + e 2A (2DE + AA'E' + E") 

d u [c; t + 2A'C fl - ie^A'G'^ - e 2A G'^ 
1 



0. 
0. 



{UH, V + Ae 2A A'H'^ + e 2 ^^) = 0. 



For (/i5) and (55) we obtain the following equations: 
d,^' + 3A'x-^'^)=0, 



e~ 2A nc^ - dg; = 0, 

4 « + 2A'ip') + 4A' X ' - e" 2A Dx + □ {E n + 2^'^') + e~ 2A UB' 



4M 3 



(B.17) 
(B.18) 

(B.19) 
(B.20) 

(B.21) 



The above equations of motion for the scalar, vector and tensor modes of perturbations are 
applicable for any gauge choice, in the main text we decide choose the longitudinal gauge 
defined by the condition B = E = = 0. 
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